Foraminifera as a model of eukaryotic genome dynamism.

In contrast to the canonical view that genomes cycle only between haploid and diploid states, many eukaryotes have dynamic genomes that change content throughout an individual's life cycle. However, the few detailed studies of microeukaryotic life cycles render our understanding of eukaryotic genome dynamism incomplete. Foraminifera (Rhizaria) are an ecologically important, yet understudied, clade of microbial eukaryotes with complex life cycles that include changes in ploidy and genome organization. Here, we apply fluorescence microscopy and image analysis techniques to over 2,800 nuclei in 110 cells to characterize the life cycle of Allogromia laticollaris strain Cold Spring Harbor (CSH), one of few cultivable foraminifera species. We show that haploidy and diploidy are brief moments in the A. laticollaris life cycle and that A. laticollaris nuclei endoreplicate up to 12,000 times the haploid genome size. We find that A. laticollaris reorganizes a highly endoreplicated nucleus into thousands of haploid genomes through a non-canonical mechanism called Zerfall, in which the nuclear envelope degrades and extrudes chromatin into the cytoplasm. Based on these findings, along with changes in nuclear architecture across the life cycle, we believe that A. laticollaris uses spatio-temporal mechanisms to delineate germline and somatic DNA within a single nucleus. The analyses here extend our understanding of the genome dynamics across the eukaryotic tree of life.IMPORTANCEIn traditional depictions of eukaryotes (i.e., cells with nuclei), life cycles alternate only between haploid and diploid phases, overlooking studies of diverse microeukaryotic lineages (e.g., amoebae, ciliates, and flagellates) that show dramatic variation in DNA content throughout their life cycles. Endoreplication of genomes enables cells to grow to large sizes and perhaps to also respond to changes in their environments. Few microeukaryotic life cycles have been studied in detail, which limits our understanding of how eukaryotes regulate and transmit their DNA across generations. Here, we use microscopy to study the life cycle of Allogromia laticollaris strain CSH, an early-diverging lineage within the Foraminifera (an ancient clade of predominantly marine amoebae). We show that DNA content changes significantly throughout their life cycle and further describe an unusual process called Zerfall, by which this species reorganizes a large nucleus with up to 12,000 genome copies into hundreds of small gametic nuclei, each with a single haploid genome. Our results are consistent with the idea that all eukaryotes demarcate germline DNA to pass on to offspring amidst more flexible somatic DNA and extend the known diversity of eukaryotic life cycles.

Forty Years of Inferential Methods in the Journals of the Society for Molecular Biology and Evolution.

We are launching a series to celebrate the 40th anniversary of the first issue of Molecular Biology and Evolution. In 2024, we will publish virtual issues containing selected papers published in the Society for Molecular Biology and Evolution journals, Molecular Biology and Evolution and Genome Biology and Evolution. Each virtual issue will be accompanied by a perspective that highlights the historic and contemporary contributions of our journals to a specific topic in molecular evolution. This perspective, the first in the series, presents an account of the broad array of methods that have been published in the Society for Molecular Biology and Evolution journals, including methods to infer phylogenies, to test hypotheses in a phylogenetic framework, and to infer population genetic processes. We also mention many of the software implementations that make methods tractable for empiricists. In short, the Society for Molecular Biology and Evolution community has much to celebrate after four decades of publishing high-quality science including numerous important inferential methods.

Somatic genome architecture and molecular evolution are decoupled in "young" linage-specific gene families in ciliates.

The evolution of lineage-specific gene families remains poorly studied across the eukaryotic tree of life, with most analyses focusing on the recent evolution of de novo genes in model species. Here we explore the origins of lineage-specific genes in ciliates, a ~1 billion year old clade of microeukaryotes that are defined by their division of somatic and germline functions into distinct nuclei. Previous analyses on conserved gene families have shown the effect of ciliates' unusual genome architecture on gene family evolution: extensive genome processing-the generation of thousands of gene-sized somatic chromosomes from canonical germline chromosomes-is associated with larger and more diverse gene families. To further study the relationship between ciliate genome architecture and gene family evolution, we analyzed lineage specific gene families from a set of 46 transcriptomes and 12 genomes representing x species from eight ciliate classes. We assess how the evolution lineage-specific gene families occurs among four groups of ciliates: extensive fragmenters with gene-size somatic chromosomes, non-extensive fragmenters with "large'' multi-gene somatic chromosomes, Heterotrichea with highly polyploid somatic genomes and Karyorelictea with 'paradiploid' somatic genomes. Our analyses demonstrate that: 1) most lineage-specific gene families are found at shallow taxonomic scales; 2) extensive genome processing (i.e., gene unscrambling) during development likely influences the size and number of young lineage-specific gene families; and 3) the influence of somatic genome architecture on molecular evolution is increasingly apparent in older gene families. Altogether, these data highlight the influences of genome architecture on the evolution of lineage-specific gene families in eukaryotes.

Testate amoebae (Arcellinida, Amoebozoa) community diversity in New England bogs and fens assessed through lineage-specific amplicon sequencing.

Testate amoebae (order Arcellinida) are abundant in freshwater ecosystems, including low pH bogs and fens. Within these environments, Arcellinida are considered top predators in microbial food webs and their tests are useful bioindicators of paleoclimatic changes and anthropogenic pollutants. Accurate species identifications and characterizations of diversity are important for studies of paleoclimate, microbial ecology, and environmental change; however, morphological species definitions mask cryptic diversity, which is a common phenomenon among microbial eukaryotes. Lineage-specific primers recently designed to target Arcellinida for amplicon sequencing successfully captured a poorly-described yet diverse fraction of the microbial eukaryotic community. Here, we leveraged the application of these newly-designed primers to survey the diversity of Arcellinida in four low-pH New England bogs and fens, investigating variation among bogs (2018) and then across seasons and habitats within two bogs (2019). Three OTUs represented 66% of Arcellinida reads obtained across all habitats surveyed. 103 additional OTUs were present in lower abundance with some OTUs detected in only one sampling location, suggesting habitat specificity. By establishing a baseline for Arcellinida diversity, we provide a foundation to monitor key taxa in habitats that are predicted to change with increasing anthropogenic pressure and rapid climate change.

Single-cell transcriptomics supports presence of cryptic species and reveals low levels of population genetic diversity in two testate amoebae morphospecies with large population sizes.

The enormous population sizes and wide biogeographical distribution of many microbial eukaryotes set the expectation of high levels of intraspecific genetic variation. However, studies investigating protist populations remain scarce, mostly due to limited 'omics data. Instead, most genetics studies of microeukaryotes have thus far relied on single loci, which can be misleading and do not easily allow for detection of recombination, a hallmark of sexual reproduction. Here, we analyze >40 genes from 72 single-cell transcriptomes from two morphospecies-Hyalosphenia papilio and Hyalosphenia elegans-of testate amoebae (Arcellinida, Amoebozoa) to assess genetic diversity in samples collected over four years from New England bogs. We confirm the existence of cryptic species based on our multilocus dataset, which provides evidence of recombination within and high levels of divergence between the cryptic species. At the same time, total levels of genetic diversity within cryptic species are low, suggesting that these abundant organisms have small effective population sizes, perhaps due to extinction and repopulation events coupled with efficient modes of dispersal. This study is one of the first to investigate population genetics in uncultivable heterotrophic protists using transcriptomics data and contributes towards understanding cryptic species of nonmodel microeukaryotes.

Updated Virophage Taxonomy and Distinction from Polinton-like Viruses.

Virophages are small dsDNA viruses that hijack the machinery of giant viruses during the co-infection of a protist (i.e., microeukaryotic) host and represent an exceptional case of "hyperparasitism" in the viral world. While only a handful of virophages have been isolated, a vast diversity of virophage-like sequences have been uncovered from diverse metagenomes. Their wide ecological distribution, idiosyncratic infection and replication strategy, ability to integrate into protist and giant virus genomes and potential role in antiviral defense have made virophages a topic of broad interest. However, one limitation for further studies is the lack of clarity regarding the nomenclature and taxonomy of this group of viruses. Specifically, virophages have been linked in the literature to other "virophage-like" mobile genetic elements and viruses, including polinton-like viruses (PLVs), but there are no formal demarcation criteria and proper nomenclature for either group, i.e., virophage or PLVs. Here, as part of the ICTV Virophage Study Group, we leverage a large set of genomes gathered from published datasets as well as newly generated protist genomes to propose delineation criteria and classification methods at multiple taxonomic ranks for virophages 'sensu stricto', i.e., genomes related to the prototype isolates Sputnik and mavirus. Based on a combination of comparative genomics and phylogenetic analyses, we show that this group of virophages forms a cohesive taxon that we propose to establish at the class level and suggest a subdivision into four orders and seven families with distinctive ecogenomic features. Finally, to illustrate how the proposed delineation criteria and classification method would be used, we apply these to two recently published datasets, which we show include both virophages and other virophage-related elements. Overall, we see this proposed classification as a necessary first step to provide a robust taxonomic framework in this area of the virosphere, which will need to be expanded in the future to cover other virophage-related viruses such as PLVs.

Foraminifera as a model of the extensive variability in genome dynamics among eukaryotes.

Knowledge of eukaryotic life cycles and associated genome dynamics stems largely from research on animals, plants, and a small number of "model" (i.e., easily cultivable) lineages. This skewed sampling results in an underappreciation of the variability among the many microeukaryotic lineages, which represent the bulk of eukaryotic biodiversity. The range of complex nuclear transformations that exists within lineages of microbial eukaryotes challenges the textbook understanding of genome and nuclear cycles. Here, we look in-depth at Foraminifera, an ancient (∼600 million-year-old) lineage widely studied as proxies in paleoceanography and environmental biomonitoring. We demonstrate that Foraminifera challenge the "rules" of life cycles developed largely from studies of plants and animals. To this end, we synthesize data on foraminiferal life cycles, focusing on extensive endoreplication within individuals (i.e., single cells), the unusual nuclear process called Zerfall, and the separation of germline and somatic function into distinct nuclei (i.e., heterokaryosis). These processes highlight complexities within lineages and expand our understanding of the dynamics of eukaryotic genomes.

Illuminating protist diversity in pitcher plants and bromeliad tanks.

Many species of plants have evolved structures called phytotelmata that store water and trap detritus and prey. These structures house diverse communities of organisms, the inquiline microbiome, that aids breakdown of litter and prey. The invertebrate and bacterial food webs in these systems are well characterized, but less is known about microbial eukaryotic community dynamics. In this study we focus on microbes in the SAR clade (Stramenopila, Alveolata, Rhizaria) inhabiting phytotelmata. Using small subunit rDNA amplicon sequencing from repeated temporal and geographic samples of wild and cultivated plants across the Northeast U.S.A., we demonstrate that communities are variable within and between host plant type. Across habitats, communities from tropical bromeliads grown in a single room of a greenhouse were nearly as heterogeneous as wild pitcher plants spread across hundreds of kilometers. At the scale of pitcher plants in a single bog, analyses of samples from three time points suggest that seasonality is a major driver of protist community structure, with variable spring communities transitioning to more homogeneous communities that resemble the surrounding habitat. Our results indicate that protist communities in phytotelmata are variable, likely due to stochastic founder events and colonization/competition dynamics, leading to tremendous heterogeneity in inquiline microeukaryotic communities.

Phylogenomic Analyses of 2,786 Genes in 158 Lineages Support a Root of the Eukaryotic Tree of Life between Opisthokonts and All Other Lineages.

Advances in phylogenomics and high-throughput sequencing have allowed the reconstruction of deep phylogenetic relationships in the evolution of eukaryotes. Yet, the root of the eukaryotic tree of life remains elusive. The most popular hypothesis in textbooks and reviews is a root between Unikonta (Opisthokonta + Amoebozoa) and Bikonta (all other eukaryotes), which emerged from analyses of a single-gene fusion. Subsequent, highly cited studies based on concatenation of genes supported this hypothesis with some variations or proposed a root within Excavata. However, concatenation of genes does not consider phylogenetically-informative events like gene duplications and losses. A recent study using gene tree parsimony (GTP) suggested the root lies between Opisthokonta and all other eukaryotes, but only including 59 taxa and 20 genes. Here we use GTP with a duplication-loss model in a gene-rich and taxon-rich dataset (i.e., 2,786 gene families from two sets of 155 and 158 diverse eukaryotic lineages) to assess the root, and we iterate each analysis 100 times to quantify tree space uncertainty. We also contrasted our results and discarded alternative hypotheses from the literature using GTP and the likelihood-based method SpeciesRax. Our estimates suggest a root between Fungi or Opisthokonta and all other eukaryotes; but based on further analysis of genome size, we propose that the root between Opisthokonta and all other eukaryotes is the most likely.

Taxon-rich transcriptomics supports higher-level phylogeny and major evolutionary trends in Foraminifera.

Foraminifera, classified in the supergroup Rhizaria, are a common and highly diverse group of mainly marine protists. Despite their evolutionary and ecological importance, only limited genomic data (one partial genome and nine transcriptomic datasets) have been published for this group. Foraminiferal molecular phylogeny is largely based on 18S rRNA gene sequence analysis. However, due to highly variable evolutionary rates of substitution in ribosomal genes plus the existence of intragenomic variation at this locus, the relationships between and within foraminiferal classes remain uncertain. We analyze transcriptomic data from 28 species, adding 19 new species to the previously published dataset, including members of the strongly under-represented class Monothalamea. A phylogenomic reconstruction of Rhizaria, rooted with alveolates and stramenopiles, based on 199 genes and 68 species supports the monophyly of Foraminifera and their sister relationship to Polycystinea. The phylogenomic tree of Foraminifera is very similar to the 18S rRNA tree, with the paraphyletic single-chambered monothalamids giving rise to the multi-chambered Tubothalamea and Globothalamea. Within the Monothalamea, our analyses confirm the monophyly of the giant, deep-sea xenophyophores that branch within clade C and indicate the basal position of monothalamous clades D and E. The multi-chambered Globothalamea are monophyletic and comprise the paraphyletic Textulariida and monophyletic Rotaliida. Our phylogenomic analyses support major evolutionary trends of Foraminifera revealed by ribosomal phylogenies and reinforce their current higher-level classification.

Genome architecture used to supplement species delineation in two cryptic marine ciliates.

The purpose of this study is to determine which taxonomic methods can elucidate clear and quantifiable differences between two cryptic ciliate species, and to test the utility of genome architecture as a new diagnostic character in the discrimination of otherwise indistinguishable taxa. Two cryptic tintinnid ciliates, Schmidingerella arcuata and Schmidingerella meunieri, are compared via traditional taxonomic characters including lorica morphometrics, ribosomal RNA (rRNA) gene barcodes and ecophysiological traits. In addition, single-cell 'omics analyses (single-cell transcriptomics and genomics) are used to elucidate and compare patterns of micronuclear genome architecture between the congeners. The results include a highly similar lorica that is larger in S. meunieri, a 0%-0.5% difference in rRNA gene barcodes, two different and nine indistinguishable growth responses among 11 prey treatments, and distinct patterns of micronuclear genomic architecture for genes detected in both ciliates. Together, these results indicate that while minor differences exist between S. arcuata and S. meunieri in common indices of taxonomic identification (i.e., lorica morphology, DNA barcode sequences and ecophysiology), differences exist in their genomic architecture, which suggests potential genetic incompatibility. Different patterns of micronuclear architecture in genes shared by both isolates also enable the design of species-specific primers, which are used in this study as unique "architectural barcodes" to demonstrate the co-occurrence of both ciliates in samples collected from a NW Atlantic estuary. These results support the utility of genomic architecture as a tool in species delineation, especially in ciliates that are cryptic or otherwise difficult to differentiate using traditional methods of identification.

Old genes in new places: A taxon-rich analysis of interdomain lateral gene transfer events.

Vertical inheritance is foundational to Darwinian evolution, but fails to explain major innovations such as the rapid spread of antibiotic resistance among bacteria and the origin of photosynthesis in eukaryotes. While lateral gene transfer (LGT) is recognized as an evolutionary force in prokaryotes, the role of LGT in eukaryotic evolution is less clear. With the exception of the transfer of genes from organelles to the nucleus, a process termed endosymbiotic gene transfer (EGT), the extent of interdomain transfer from prokaryotes to eukaryotes is highly debated. A common critique of studies of interdomain LGT is the reliance on the topology of single-gene trees that attempt to estimate more than one billion years of evolution. We take a more conservative approach by identifying cases in which a single clade of eukaryotes is found in an otherwise prokaryotic gene tree (i.e. exclusive presence). Starting with a taxon-rich dataset of over 13,600 gene families and passing data through several rounds of curation, we identify and categorize the function of 306 interdomain LGT events into diverse eukaryotes, including 189 putative EGTs, 52 LGTs into Opisthokonta (i.e. animals, fungi and their microbial relatives), and 42 LGTs nearly exclusive to anaerobic eukaryotes. To assess differential gene loss as an explanation for exclusive presence, we compare branch lengths within each LGT tree to a set of vertically-inherited genes subsampled to mimic gene loss (i.e. with the same taxonomic sampling) and consistently find shorter relative distance between eukaryotes and prokaryotes in LGT trees, a pattern inconsistent with gene loss. Our methods provide a framework for future studies of interdomain LGT and move the field closer to an understanding of how best to model the evolutionary history of eukaryotes.

Newly designed foraminifera primers identify habitat-specific lineages through metabarcoding analyses.

Foraminifera include diverse shell-building lineages found in a wide array of aquatic habitats from the deep-sea to intertidal zones to brackish and freshwater ecosystems. Recent estimates of morphological and molecular foraminifera diversity have increased the knowledge of foraminiferal diversity, which is critical as these lineages are used as bioindicators of past and present environmental perturbation. However, a comparative analysis of foraminiferal biodiversity between their major habitats (freshwater, brackish, intertidal, and marine) is underexplored, particularly using molecular tools. Here, we present a metabarcoding survey of foraminiferal diversity across different ecosystems using newly designed foraminifera-specific primers that target the hypervariable regions of the foraminifera SSU-rRNA gene (~250-300 bp long). We tested these primer sets on four foraminifera species and then across several environments: the intertidal zone, coastal ecosystems, and freshwater vernal pools. We retrieved 655 operational taxonomic units (OTUs); the majority of which are undetermined taxa that have no closely matching sequences in the reference database. Furthermore, we identified 163 OTUs with distinct habitat preferences. Most of the observed OTUs belonged to lineages of single-chambered foraminifera, including poorly explored freshwater foraminifera which encompass a clade of Reticulomyxa-like forms. Our pilot study provides the community with an additional set of newly designed and taxon-specific primers to elucidate foraminiferal diversity across different habitats.

Macronuclear development in ciliates, with a focus on nuclear architecture.

Ciliates are defined by the presence of dimorphic nuclei as they have both a somatic macronucleus and germline micronucleus within each individual cell. The size and structure of both germline micronuclei and somatic macronuclei vary tremendously among ciliates. Except just after conjugation (i.e. the nuclear exchange in their life cycle), the germline micronucleus is transcriptionally inactive and contains canonical chromosomes that will be inherited between generations. In contrast, the transcriptionally active macronucleus contains chromosomes that vary in size in different classes of ciliates, with some lineages having extensively fragmented gene-sized somatic chromosomes while others contain longer multigene chromosomes. Here, we describe the variation in somatic macronuclear architecture in lineages sampled across the ciliate tree of life, specifically focusing on lineages with extensively fragmented chromosomes (e.g. the classes Phyllopharyngea and Spirotrichea). Further, we synthesize information from the literature on the development of ciliate macronuclei, focusing on changes in nuclear architecture throughout life cycles. These data highlight the tremendous diversity among ciliate nuclear cycles, extend our understanding of patterns of genome evolution, and provide insight into different germline and somatic nuclear features (e.g. nuclear structure and development) among eukaryotes.

Epigenetic influences of mobile genetic elements on ciliate genome architecture and evolution.

Mobile genetic elements (MGEs) are transient genetic material that can move either within a single organism's genome or between individuals or species. While historically considered "junk" DNA (i.e., deleterious or at best neutral), more recent studies reveal the potential adaptive advantages MGEs provide in lineages across the tree of life. Ciliates, a group of single-celled microbial eukaryotes characterized by nuclear dimorphism, exemplify how epigenetic influences from MGEs shape genome architecture and patterns of molecular evolution. Ciliate nuclear dimorphism may have evolved as a response to transposon invasion and ciliates have since co-opted transposons to carry out programmed DNA deletion. Another example of the effect of MGEs is in providing mechanisms for lateral gene transfer (LGT) from bacteria, which introduces genetic diversity and, in several cases, may drive ecological specialization in ciliates. As a third example, the integration of viral DNA, likely through transduction, provides new genetic materials and can change the way host cells defend themselves against other viral pathogens. We argue that the acquisition of MGEs through non-Mendelian patterns of inheritance, coupled with their effects on ciliate genome architecture and persistence throughout evolutionary history, exemplify how the transmission of mobile elements should be considered a mechanism of transgenerational epigenetic inheritance.

Examining the Relationship Between the Testate Amoeba Hyalosphenia papilio (Arcellinida, Amoebozoa) and its Associated Intracellular Microalgae Using Molecular and Microscopic Methods.

Symbiotic relationships between heterotrophic and phototrophic partners are common in microbial eukaryotes. Among Arcellinida (Amoebozoa) several species are associated with microalgae of the genus Chlorella (Archaeplastida). So far, these symbioses were assumed to be stable and mutualistic, yet details of the interactions are limited. Here, we analyzed 22 single-cell transcriptomes and 36 partially-sequenced genomes of the Arcellinida morphospecies Hyalosphenia papilio, which contains Chlorella algae, to shed light on the amoeba-algae association. By characterizing the genetic diversity of associated Chlorella, we detected two distinct clades that can be linked to host genetic diversity, yet at the same time show a biogeographic signal across sampling sites. Fluorescence and transmission electron microscopy showed the presence of intact algae cells within the amoeba cell. Yet analysis of transcriptome data suggested that the algal nuclei are inactive, implying that instead of a stable, mutualistic relationship, the algae may be temporarily exploited for photosynthetic activity before being digested. Differences in gene expression of H. papilio and Hyalosphenia elegans demonstrated increased expression of genes related to oxidative stress. Together, our analyses increase knowledge of this host-symbiont association and reveal 1) higher diversity of associated algae than previously characterized, 2) a transient association between H. papilio and Chlorella with unclear benefits for the algae, 3) algal-induced gene expression changes in the host.

DAPI staining and DNA content estimation of nuclei in uncultivable microbial eukaryotes (Arcellinida and Ciliates).

Though representing a major component of eukaryotic biodiversity, many microbial eukaryotes remain poorly studied, including the focus of the present work, testate amoebae of the order Arcellinida (Amoebozoa) and non-model lineages of ciliates (Alveolata). In particular, knowledge of genome structures and changes in genome content over the often-complex life cycles of these lineages remains enigmatic. However, the limited available knowledge suggests that microbial eukaryotes have the potential to challenge our textbook views on eukaryotic genomes and genome evolution. In this study, we developed protocols for DAPI (4',6-diamidino-2-phenylindole) staining of Arcellinida nuclei and adapted protocols for ciliates. In addition, image analysis software was used to estimate the DNA content in the nuclei of Arcellinida and ciliates, and the measurements of target organisms were compared to those  of well-known model organisms.The results demonstrate that the methods we have developed for nuclear staining in these lineages are effective and can be applied to other microbial eukaryotic groups by adjusting certain stages in the protocols.